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Abstract 


We retrace the recent history of the Umbral Calculus. After studying the classic 
results concerning polynomial sequences of binomial type, we generalize to a 
certain type of logarithmic series. Finally, we demonstrate numerous typical 
examples of our theory. 

Une Presentation des Contributions Recents au Calcul des Diffe¬ 
rences Finies 

Nous passons en revue ici les resultats recents du calcul ombral. Nous 
nous interessons tout d’abord aux resultats classique appliques aux suites de 
polynomes de type binomial, pius elargions le champ d’etude aux series loga- 
rithmiques. Enhn nous donnons de nombreaux exemples types d’application de 
cette tlieorie. 
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1 Polynomials of Binomial Type 


Much of the calculus of finite differences leans upon certain analogies between 
two linear operators defined on functions f(x) ranging over a field K of charac¬ 
teristic zero, and taking values in K: 


1. The derivative operator D, defined as 


D f{x) = f'{x), 


and 


2. The forward difference operator A, defined as 


A/(x) = f{x + 1) - f(x). 


1.1 The Powers of x 

The sequence of polynomials 1 ] p n (x) = x n is related to the derivative operator 
by the elementary property 


Dp„(x) = np n -i{x). (1) 

Such sequences of polynomials have been called Appell. The sequence x n is 
distinguished among all Appell sequences by the further property 

Pn ( 0 ) = &n0 


for all n > 0. 

1 When one speaks of a sequence of polynomials, one generally means a sequence p n (x) 
indexed by a nonnegative integer n such that the degree of each polynomial p n {x) is exactly 
equal to n. 
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Furthermore, p n {x) = x n plays a key role in Taylor’s theorem 

„ , ^ ^ "/(*) ^a n D n f(x) 

f(x + a) = ^ Pn{a) , = 2^- - 

n —0 * n —0 

for all a £ K and for all suitable functions f(x). 

Finally, we note that the sequence p n (x ) is of binomial type. That is to say, 
it satisifies the binomial theorem 


p n (x + a) = ^2 (j^jPk(a)Pn-k(x) 


( 2 ) 


for all a £ K and all n > 0. 


1.2 The Lower Factorial 


Define another sequence of polynomials called the lower factorial powers as 
follows 

Po{x) = (z)o = 1 


and 


Pn{x) = (x)n = x(x - 1) {x — 2) ■ • ■ (x — U + 1) 

for n > 0. It can be interpreted as the number of injections from an n-element 
set to an x-element set.Q 


The lower factorial is related to the forward difference operator by the prop¬ 
erty 

A p n {x) = np n - i(x). 

Such sequences are called Shejfer sequences for the forward difference operator. 
The lower factorial is distinguished among all Sheffer sequences by the further 
property 

Pn (b) — 3n0 

for n > 0. 


Furthermore, for suitable functions f(x) and all constants a £ K, we have 
Newton’s formula 


OO 

f(x + a) = ^2(a) n 

n =0 


A n f(x) 
n\ 


Lastly, Vandermonde’s identity states that the sequence p n {x) = (x) n is of 
binomial type: 

OO 

(x + a)n = 

k =0 


(a)k(x)n- k . 
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1.3 The General Case 


Our initial problem will be that of carrying this analogy as far as reasonably 
possible in order to make the parallel of D and A into an accidental special case. 

To this end, we begin by classifying all sequences of polynomials of binomial 
type, that is, all sequences of polynomials p n (x) (with deg (p n (x)) = n satisfying 
eq-l 

OO 

Pn{x + a) = 

k—0 


Pk(a)p n -k(x). 


Proposition 1 A sequence of polynomials p n {x) of binomial type satisfies 

Pn ( 0 ) = bnO 


for all n > 0. 


Proof: By eq. |] with n = a = 0, po(0) = 1. Now, for n positive, 

0 = p n (x) - p n (x) = ^ (j^jPk(Q)Pn-k{x). 

However, pk{x) is a basis for the vector space of polynomials, and in character¬ 
istic zero the binomial coefficients indicated are never zero, so Pfc(0) = 0.D 

Now, to each sequence of binomial type p n {x), we associate a linear operator 
Q defined as 

Qp n {x) = np n -i(x). (3) 

For the moment, the operator Q is defined only for polynomials; later (Sec¬ 
tion ^), we shall extend its domain to all more general types of functions. 

By iteration, we obtain 


Q k Pn(x ) = ( n) k Pn-k{x ). 

Recalling that ((() = ( n)k/k\ , we can write the binomial identity (eq. |^) in the 
form 

Pj (x + a) = £><«) A|M 

k =0 

for j > 0. 
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Now, let p(x) be any polynomial. Given its expansion p(x) = ^°1 0 c jPj( x ) 
in terms of the basis pj ( x ), we have 

OO 

p(x + a) = CjPj (x + a) 

k= 0 

OO / OO \ 

= ^2pk(a)Q k Y/zpAx) 
k—0 \j~ 0 / 

OO 

= Y pk ^~k[ p ^- ( 4 ) 

k =0 

This calculation for all polynomials p{x) can be recast into a more elegant form 
by introducing the shift operator 

E a f{ x ) = f(x + a). 

We then can write eq. [| as the operator identity 

OO 

E a = J2 p k( a )Q k / k '- ( 5 ) 

k= 0 

which we shall call the Taylor's Formula for the sequence p n (x). Hence, New¬ 
ton’s formula is nothing more than Taylor’s formula for the lower factorials. 
Again, eq. has been proven only for operations on polynomials for the time 
being; however, we shall see (Theorem ^|) that its domain of validity is consid¬ 
erably greater. 

The right side of eq. [ 5 ] obviously commutes with Q. Thus, 

QE a = E a Q 

for all constants a. Any linear operator T with the property that 

TE a = E a T (6) 

for all scalars a will be said to be a shift-invariant operator. A linear shift- 
invariant operator is called a delta operator if its kernel is exactly the field of all 
constants K. 

We have thus proved all of Theorem [[] and the hard part of Theorem |[ 

Theorem 1 Let p n {x) he a sequence of polynomials of binomial type. Then 
there exists a unique delta operator such that 

Qp n {x) = np n -i(x) 

for n > (].□ 
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Next, we have the converse. 


Theorem 2 Let Q be a delta operator. Then there exists a unique sequence of 
polynomials of binomial type p n (x) satisfying 

Qp n {x) = np n - i(x) 


for n > 0. 


Proof: We can recursively define p n {x ) subject to the requirements that 


• Po(x) = 1 , 

• For n positive, p„fx ) is in the inverse image of Q on np n _i(x), and 

• For n positive, p„(0) = 0. 

It remains only to show that the resulting sequence of polynomials obeys eq. || 
By the reasoning before Theorem [i], we again have eq. [|. Finally, setting p{x) = 
Pn(x) yields the binomial theorem.□ 

Thus, we see that all of the features of the analogy between the operators D 
and A, and the polynomial sequences x n and ( x) n are shared by all other delta 
operators, and their corresponding sequences of binomial type. 

Sequences of polynomials of binomial type are of frequent occurence in com¬ 
binatorics^] probability, statistics, function theory and representation theory. 


1.4 The Abel Polynomials 

We note the following striking example of a sequence of polynomials of binomial 
type. Given a constant b , define Abel’s polynomials [j] as follows 

A 0 {x ; b) = 1 


and for n positive 

A n (x; b) = x(x — nb) n ~~ l . 

2 In combinatorics, sequences of binomial type enumerate the number of functions from an 
n-element set to an ai-element set enriched with some mathematical structure on each fiber. 
For details see |H- 

3 For example, the Abel polynomial A n (x\— 1) counts the number of reluctant functions 
from an n-element set to an x-element set. A reluctant function is a function enriched with a 
rooted labelled forest on each fiber. 


6 



It is easily verified that, setting Q = E b D = D E b , we have QA n (x\b) = 
nA n -i(x;b). Hence, Theorem implies at once Abel’s identity 

OO 

(x + a)(x + a — nb) n_1 = ^ 

fc=o 

Furthermore, Taylor’s formula for the sequence of Abel polynomials gives 

p(x + a) = £ a(a ~ fc f )fc " 1 D k p(x + kb). 
k—0 

For example, for p(x) = x n , we obtain 

(x + a) n =jr ° ( °~^ \ n) k (x + kb) n ~ k 

k -0 

which is not a priori obvious. 

In order to devise more examples, it will be helpful to first study shift- 
invariant operators in greater detail. 




2 Shift-Invariant Operators 


For convenience, we introduce the following bit of notation. For any function 
f{x), we write 

(/0))(o) = /(°) = [f( x ) J*=o- (7) 

The linear functional ()^ 0 ^ is called the augmentation or evaluation at zero. 

Let p n {x) be a sequence of binomial type and let Q be its delta operator. 
Taylor’s theorem for p n (x) can be written in terms of the augmentation (inter¬ 
changing the roles of x and a) 


E a p(x) = £>*(*) 

k=0 


(E a Q k p(x)) (0) 
k\ 


Now, let T be any shift-invariant operator, and p{x) any polynomial. We have 


E a Tp{x) 


TE a p{x) 

OO 

J2T Pk (x) < E a Q k p(x)) {0] /k\. 

k—0 
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Again exchanging the role of x and a, therefore 


Tp(x + a) 


OO 


E 




Setting a = 0, we obtain the operator identity 


T = £c fc Q fc /*! 

k=0 


( 8 ) 


where c k = (Tp k (x))^ 0 y 
Thus, we have proven 

Theorem 3 (Expansion Theorem) Given a shift-invariant operator T and 
a delta operator Q associated to a sequence of binomial type p n (x ) there is a 
unique expansion given by eq. |^.n 


Actually, some mild continuity conditions must be imposed, but we will refer the 
reader to the bibliography for such technical (and ultimately trivial) questions 
of topology. 

The expansion theorem is to operators essentially what the generalized Tay¬ 
lor’s theorem is to polynomials. Together they allow the expansion of either 
operators or polynomials any basis. 

The above result can be restated in more elightening terms. Let K[[t]\ be 
the ring of all formal series (sometimes called Hurwitz series ) in the variable t 
with coefficients in K 

f(t) = Co + c\t + c 2 t 2 /2\ + c 3 f 3 /3! H-. 


Then we have 

Theorem 4 (Isomorphism Theorem) 1. For every delta operator Q, one 
has the following isomorphism of the ring of formal power series K [[£]] 
onto the ring of shift-invariant operators 

OO 

/(*) ~ 5>*Q fc /*! =/(Q)- (9) 

k=0 


2. Furthermore, let (f> be an isomorphism of the ring of formal power series 
with the ring of shift-invariant operators. Then there exists a delta oper¬ 
ator Q such that the isomorphism (j> is of the form eq. [ij. 
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Proof: Part one is immediate from the Expansion Theorem. Part two is 
classical result concerning formal series where Q = </>(£).□ 

In particular for Q = D, we notice that all shift-invariant operators are 
formal power series in the derivative. That is, the ring of shift-invariant opera¬ 
tors is if[[D]]. Moreover, delta operators T = T(D) are formal power series in 
the derivative inwhich the coefficient of D° is zero, and the coefficient of D 1 is 
nonzero. Such series are known as delta series. 

For example, the shift operator is given by the formal power series 


OO 


exp(aD) = ^ a n D n /n\. 


E a 


As another example, let K = C be the complex field, and define the Laguerre 
operator by 


The Lagerre operator is obviously a shift-invariant operator since 



( 10 ) 



Moreover, 




'OO 


I e v n\dy 
o 


and hence by the Expansion Theorem (Theorem ||) , we infer that L is a delta 
operator and that 


L = -D - D 2 - D 3 
= D/(D-7) 


( 11 ) 
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where I = D° is the identity operator. We will later (Section |4.2| ) compute the 
sequence of binomial type associated with the Laguerre operator. 


As another example, let W be the Weierstrass operator defined as 

1 r°° 

Wf(x) = —= / f{y) exp(-(a; - yf/2)dy. 

V J — oo 


Again, we calculate c n = (.Lx n ), 0 y First, note that for n odd, the integrand in 
question is an odd function, so the integral from —oo to oo is zero. For n = 0, 
it is easier to compute Cq. 


c n = 


^ /»oo pOO 

— / exp(—(a; 2 + y 2 )/2)dxdy 

i p2ir poo 

— / exp(— r 2 /2)r dr d6 
27T J o J o 

/»oo 

/ exp(— r 2 /2)rdr 

Jo 


exp (—u)du 


= 1. 


In general, by integration by parts, we have the following recurrence for c n 
1 r°° 

Cn = —f= / y n enp(-y 2 /2)dy 

V J —oo 

= -±= (- [exp(—y 2 /2)y n ” 1 ] ^ + (n - 1) J°° y n ~ 2 exp(-r/ 2 /2)dy) 
= 0 + (n — l)c„_ 2 - 

Thus, for n even, 


so that 


c n = [n — 1) (n — 3) • • • 5 • 3 • 1 = 


(n/2)\2 r ‘ 


n=0 


D 2 ™ 
n\2 n 


exp(D 2 /2). 


From the Isomorphism Theorem (Theorem Q) , it follows that a shift-invariant 
operator has a unique inverse if and only if T\ ^ 0. For example, the Weierstrass 
operator W has a unique inverse. 
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Lastly, we compute the Bernoulli operator 


Jf(x) 


rx+l 


f(y)dy 


where we have 


Thus, 


(Jx n ) i0 ) ~ 



,.n+1 


1 


n - 


y =o 


1 

n + 1 


OO 


k—0 


D k 

(fc + 1)! 


e D — / A 
D = D' 


We close this section with a result which may be considered to be funda¬ 
mental. 


Theorem 5 (Fundamental Theorem of the Calculus of Finite Differ¬ 
ences) Let p n {x) be the sequence of binomial type associated to a delta operator 
Q(D). LetQ^ l \t) be the inverse formal power series ofQ(t). That is, suppose 
that <3(<3 (-1) (^)) = ^(Q{t)) = t. Then the exponential generating function 

for p n {x) is 

OO 

Y,Pn{x)t”/n\ = exp(xQ (-1) (f)). 

n—0 


Since this result is so important, we give two proofs. A third proof found in 
|5) is more powerful and applies equally to symmetric functions. 

Proof 1: By Proposition [I] and eq. ||, the left hand side is characterized by 
the property 

(Q(D) k p{x,t)) {Q) =t k , 

so it will suffice to verify the right hand side obeys this property equally. Let 
q n k denote the coefficients of Q(D) fe = 9 nfcD n /n!. Then 


(<2(D) fe exp(a;Q ( 1 ] (t))\ 

\ / ( 0 ) 


= ^<Q(D) fc ; r”) (0) (Q(- 1 )(t))"/n! 

n =0 
oo 

= E dnk{.Q [ ~ l) {t)) n /n\ 
n —0 

= t k .a 
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Proof 2: By Taylor’s Formula (eq. ||), 

OO 

exp(oD) = 5 ^Pk(a)Q k /k\. 
k=0 

Substituting x for a, and for D, we have the desired result.□ 


3 Extension to Logarithmic Series 


3.1 The Harmonic Logarithms 


We shall now extend the domain of every shift-invariant operator to a more 
general domain of formal series originating from the Hardy fieldQ C called the 
Logarithmic Algebra consisting of all expansions of real functions in a neigh¬ 
borhood of infinity in terms of the monomials 

t a = x a ° log(cc) ai log(log(a;)) a2 • • • 

for all vectors of integers^ where only a finite^] number of the on are different 
from zero. 


That is C is equivalent to the set of all formal sums 

a. 

where the sum is over vectors of integers a with finite support, and for all 
nonnegative integers n, all integers ao, aq,..., a„_i, and all integers f3, there 
exists finitely many a n > j3 such that there exists integers a n + 1 , a n + 2 , • ■ ■ (only 
finitely many different from zero) such that b a ^ 0. 

Another characterization of C |(], Theorem 5.4A.4] is that C contains a dense 
subset which is the smallest proper field extensions of the ring of polynomials 
with real coefficients such that the derivative D is a derivation of the entire field, 
and the antiderivative D _1 is well defined up to a constant of integration. 


It will be noted that the derivative of such an expression is awkward 


Dr = a 0 £ 


l (aO- 1 ’ Q D Q 2> Q 3"0 


OL\i 


.(ao-l.ai-l,. 


* 3 ") 


Oi^l 


(ot 0 -!,<*!-1,0(2-1,013...) 


4 No knowledge of Hardy fields nor of asymptotic expansions is expected of the reader. 

5 We could achieve far greater generality by allowing a to be a vector of reals as in |bj. The 
definitions in p] provide one with the basic tools to carry out these calculations. 

6 Since all the vectors we deal with here have finite support, we will usually adopt the 
convention of not writing the infinite sequence of zeroes they all end with. Thus, we write (1) 
for the vector (1,0,0,0,... 
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Moreover, the formula for the antiderivative of such a function is not known 
in general, and even in the case ao 7 ^ 1 when it is known, it is given by a 
horrendous expression || p. 87]. Since the monomials £“ are so unwieldy and 
do not even begin to compare to the powers of x in terms of ease and simplicity 
of calculation, we are forced to resort to another basis for the vector space C 
which as will be seen is the true logarithmic analog of the sequence x n . This 
basis is called the harmonic logarithm A“(*). 


For n an integer, and a = (cti, 02 ,...) a vector of integers with only a finite 
number of nonzero entries, the harmonic logarithm of degree n and order a is 
defined by the series 


An(z) =Y1 

P 


s{~n, pi) 
L-n]! 


OO 

jPJ ^Pi—Pi+i Q-i 1; • • • 5 CXi Pi T 1) 

_*=1 


g(n),(a-p) 


where 

1. The linear partition p = (pi,p 2 , ■ ■ •) is a weakly decreasing sequences of 
nonnegative integers, 

2. The Roman factorial^ rn\ is given by the formula 

lU j n\ for n a nonnegative integer, and 

L ’ [ (—l ) rl_1 /(—n — 1)! for n a negative integer, 

3. For any integer n negative, zero or positive, and for any nonnegative in¬ 
teger k , the Stirling number s{n 1 k) is the coefficient of y k in the Taylor 
series expansion of (y) n = r(y + 1)/T(y — n + 1), and 

4. The elementally symmetric function e n (xi,X 2 , • • ■) is defined the generat¬ 
ing function 

OO OO 

y, e n (xi,X 2 , ■ ■ ■)y n = J|(l + x k y). 

71=0 fc= 1 


It may be helpful to point out several important special cases. 


1. The harmonic logarithms of degree zero are given by 

Ag = ^(°)>“ = (logar)“(logloga:)2 • • ■ • 


2. The harmonic logarithm of order a = (— 1) and degree one is the 

logarithmic integral li(a;) = f 0 dt/log t. 

'After Steve Roman. |l5|, |l6[ 
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3. The harmonic logarithms of order a = (0) are given by 


Ai 0) (*) 


x n for n > 0, and 
0 for n < 0. 


4. 


5. 


(log x — 1 — 2 — ''' — n) for n > 0, and 

for n < 0. 


The harmonic logarithms of order a = (1) are given by 
X ( °\x) = 

The harmonic logarithms of order a = (2) are given by 
A^ 0) (x) = 


X n (logx) 2 — X n (j + I +-b —) 

~^x n ^ (1 + |) Tj’ ■ • -I- " for n > 0, and 

n (-*■ + 2^ ' + n)] 


2x n 


l°g( a; ) “ 1 - 3- 


for n < 0. 


6. The harmonic logarithms of order a = ( t ) where t is a nonnegative integer^ 
are given by the following formulas 





[n] \x n+z T{z + l)/T{z + n + 1) 
x n ln]\(xD)_ n (\ogx) t 

t 

x n E [n]\(t)ks(-n,k)(\ogxY~ k . 
fe=o 


Now we can give the association between the derivative and the harmonic 
logarithms; the harmonic logarithms behave under derivation exactly like the 
powers of x once the ordinary factorial n! is replaced by the Roman factorial 

Mi- 


Theorem 6 For all integers n, and vectors of integers with finite support a, 
we have 


D K+i(x) = [n + 1] A“(ir) 


( 12 ) 


where 


, Ini! f n forn^O, and 

W = ISVTj! = { 1 for n = 0. 
More generally, for any nonnegative integer k, 


D k xa _ + 

U A n+k ~ ^-| j A n\ X >- 

8 Our original paper ]l2| dealt exclusively with this case. 
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Proof: By induction, it will suffice to demonstrate eq. |l2]. The following 
proof is a direct application of the recursion for the Stirling numbers. (See 

0 -) 

DA “+iO) 


= E 


= E 


s{—n — l,pi) 
L ~n - 11! 

s(-n- l,pi) 

L-«- 11 1 


f'Pi — pi+l (^hi 1; * * * ) Clij Pi “b 1) 

. 1=1 
’ oo 

€pi—pi -\-1 (^25 1? • • • j Pi 1) 


Y)£(n+l),{cx-p) 


. 2—1 

OO 


(a + 1 - p/ , 


fe=i 


= E 


(a + l)s(—n - 1, pi) + s(—n - 1, pi) 

l-n~ 11 ! 

I-’ 

’ OO 

| & pi— Pi+\ {pii 1} • * ■ ) Pi T 1) 


. 1=1 


g( n ),( a ~P) 


y- s (~ n i Pi) 

p L 1 

= L n + 11 A n(^)- D 


^ | & Pi— Pi+\ {p-ii Ch 1) • • • ) Pi "b 1) 


i=l 


^(n.),(a-p) 


Hence, if we denote by C + the closure of the span of the harmonic logarithm 
A“(x) with a 7 b 0, then we notice that D restricted to C + is a bijection. In 
particular, the antiderivative is given by 

D ^K( X ) = A“ + 1 (x)/ [n + 1] . 


We can thus characterize the harmonic logarithms by the relationship. 
A“(x) = [n}\ D" n £(°l“ = [nl!D- n A((. 


Now, let us dissect the logarithmic algebra even more finely. Let CP denote 
the logarithmic series involving A“(x) for all integers n and a fixed a. Eviden- 
tally, C + is the direct sum of the C a for a 7 ! (0). 


Theorem 7 (Roman Modules) All of the vector spaces C a for a 7 ^ (0) are 
naturally isomorphic as differential vector spaces. 
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Proof: The isomorphism is given by the skip operator 
skip a/3 A^(:r) = A“(x).D 


In view of the proceeding theory, most of our calculations can be made in 
Zb 1 ) w hich is called the Roman module after Steve Roman. The Roman module 
consists of series of the form 


f(x) = 


-1 

= E 

j=-oo 


CjX 3 + 


E C 3® J 

3=0 


log x — 1 - 

2 j 


= E 




(i) 


0 ); 


] = ~oo 


that is, there are a finite number of terms in which logx appears, and an infinite 
sereies in inverse powers of x. Thus, the Roman module is fairly convenient for 
our calculations. All of our concrete examples will be drawn from the Roman 
module. 


However, /b°) is not isomorphic to the Roman module; it is isomorphic to 
the algebra of polynomials C[x]. 


3.2 Shift-Invariant Operators 

A linear operator Q will be called shift-invariant if not only does it commute 
with the shift operator E a = Yln=o «"D n /n\ for all a, but also that it commutes 
with the skip operators skip aj g from the logarithmic algebra C to C a 

skip a/3 A^(x) = <5/3 7 A“(x) 

for (3 ft (0). 

Any shift-invariant operator on C is obviously shift-invariant (in the sense of 
eq. |) when restricted to /b°) = C[x]. Conversely, any shift-invariant operator 
Q on C[x] can be expanded as a formal series in the derivative Q(D), and thus Q 
can be extended to all C. However, are these the only shift-invariant operators 
which act on £? On C + 1 

Theorem 8 (Characterization of Shift-Invariant Operators) 

(1) The algebra of shift-invariant operators on C + is naturally isomorphic 
to the algebra of Laurent serfe^] in the derivative C (D). 

(2) The algebra of shift-invariant operators on C is naturally isomorphic to 
the algebra of formal power series in the derivative C[[D]]. 

9 A Laurent series is a formal series of the form 

OO 

f(x) = ^ CiX 1 

i=n 

where n may be any constant positive or negative. 
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Proof: It will suffice to prove part one, since the only Laurent series in the 
derivative which are well defined on are those which involve no negative 
powers of the derivative; that is, members of the algebra of formal power series 
in the derivative. 


Since the derivative commutes with skip a ^ and E a , all Laurent series in the 
derivative are shift-invariant operators on C + . 

Conversely, suppose that Q is a shift-invariant operator. Q is determined by 
its actions on the harmonic logarithms of order a for any particular a 7 ^ ( 0 ). 

n 

QK{X)= £ CnkK{x). 

k =—00 


Next, since Q commutes with E a for all a, it also commutes with D n for all n. 
Thus, 

* LH! 




j=-oo 


QV n \%(x) 

I h 1 ! k ~ n 

[k-n 1! Cn ~ k ’ jX ^ x) - 


J=—oo 


Equating coefficients of A“(x), and setting a = b, we have 


c 0 j — 


J 

k 


Ck,j-\-k 


where the Roman coefficient 


j 

k 


is defined to be |_j]!/ [k~\ ! [j — k ]!. 


Hence, Q is determined by the coj, and therefore equals the Laurent series 
£™_oo cojD J / |_j~|! where m is the largest j such that Coj ^ 0 .D 

Since the algebra of Laurent series is a field, we immediately derive that all 
shift-invariant operators on C + are invertible. Thus, 


Corollary 1 (Differential Equations) (1) All linear differential equations 
with constant coefficients on C + have a unique solution. That is, any equa¬ 
tion of the form 

f(D)p(x) = q{x) 

(where q(x) £ C + ) has a unique solution p(x) £ C + . 

(2) All linear differential equations on C have solutions. Of these, one par¬ 
ticular solution can be naturally chosen as the canonical one. 

Proof: We have just proven 1. The proof for 2 (see j?], §2.4.4A]) relies on the 
use of the projection maps skip^/TD 
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3.3 Binomial Theorem 


Any shift-invariant operator acts on each level C a in the same way. For example, 
consider the logarithmic version of the binomial theorem. 


E a \«{x) = Y J a k V k K{x)/k\ = Y J 

k=0 k—0 


n 

k 


a k \°_ k (x). 


(13) 


Several special cases are of particular interest. For a = (0), we have the standard 
binomial theorem. For a = (1), n positive, and x = 1, we have i P- 97] 


(1 + a) n log(l + a) 


/ 1 

i\ 

( 1 

1 \ 

i + - + •• 

■ • H— ) 

— na I 1 + - + • ■ 

* + ) 

V 2 

n 

V 2 

n — 1 / 


~\~d 


n 

2 

n+1 


a ^ l + \ + -" + ^2 


2 \n — 2 


_n—2 _ n—1 

a — na 


n 


+ a 


n+2 


n 

-2 


+ ■ 


k -1 

Similarly, for a = (0,0,..., 0,1) and n = 0, 

n 

log log • ■ ■ log(x + a) 

n 

= loglog • • ■ log(x) + ( —l)Pi + -+Pfc- 2 + 2 Pfc-i + l 

Pl>P2>->Pfe=l 


X 


k—2 

U e Pj~Pj +1 (1) ■ ■ • ) Pj — 1) 
1=1 


53 zU ( n ^ 

n =0 


This framework immediately gives rise to logarithmic versions of Taylor’s 
theorem and the expansion theorem in terms of the logarithmic generalization 
of evaluation at zero or augmentation. The augmentation^ of order a is defined 
to be the linear functional () from the logarithmic algebra C to the complex 
number C such that 

( A n(x)) a = S a/ 3 S n0 . 


10 The restriction of the augmentation of order (0) to is the evaluation at zero mentioned 

in eq. [71 The other augmentations can be profitably be rewritten in terms of the residue at 
zero. 
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Theorem 9 (Logarithmic Taylor’s Theorem) Let p{x) £ C + . Then we 
have the following expansion of p(x) in terms of the harmonic logarithms. 


p(x) 


E E 

a^(0) n=oo 


(D n p(x)) a 

N! 


K( x )- 


( 14 ) 


Proof: It suffices to show that eq. 0 
Theorem || states that eq. 0 holds for the 


holds for a basis of C + . However, 
basis of harmonic logarithms.□ 


Theorem 10 (Logarithmic Expansion Theorem) (1) Let Q(D) he a shift- 
invariant operator on C + . Then we have the following expansion of Q(D) in 
terms of the powers of the derivative 


OO 

Q(D)= J2 

n=—oo 


(Q(D) K{x)) a 

M! 


where a / (0). 

(2) Let Q(D) be a shift-invariant operator on C. Then we have the following 
expansion of Q(D) in terms of the powers of the derivative 

n —0 


where a is any vector of integers with finite support. 


Proof: Theorem is immediately for Q(D) = D n from Theorem !□ 

We now see that the sequence of harmonic logariths A“ ( x ) is in some sense 
(to be made clear later) of “binomial type,” and that in this sense the harmonic 
logarithms are associated with the derivative operator. 


3.4 The Logarithmic Lower Factorial 


We can now extend to the logarithmic domain the analogy between the deriva¬ 
tive operator D and the forward difference operator A which has been noted for 
polynomials. 

Since £^0) is the algebra of polynomials, one would rightly expect that the 
logarithmic version Pn(x) = (x)“ of the lower factorial function p n (x ) = ( x) n 
would include the subsequence 


Pn\x) = (x)i 0) = 


(x) n = x(x — 1) • • • (x — n + 1) for n nonnegative, and 
0 for n negative. 
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It has been known for a long time that, setting for n negative, 


Pn{x) = ( x)n = l/(x+l)(x + 2 ) • ■ ■ (x - Tl), 


one has 

A p n (x) = n(x) n - 1 - 

Given this preliminary information]^] we guess that when n is negative, l/(x + 
1 )■■■ (x — n) is to the forward difference operator as x n is to the derivative. 
Thus, we write 


Pn ] (x) = (x)^l) = 1/(x + 1)(x + 2) ■ ■ ■ (x + n) 


for n negative. 

This leave open the question of what is the forward difference version of the 
logarithm Aq 1 '^) = log a; and of the harmonic logarithms 

W(x)=x n (logz-l-i- 

for n positive. In other words, what is the formal solution p^ ( x ) = (x)q 1 ^ of the 
difference equation A p(x) = l/(x + 1)? This question was first formulated by 
Gauss; in the present context, it can be dealt with very easily. By Corollary [l| 
all such difference equations have a unique solutions. 

In fact, we can easily calculate (x)q 1 '*. Since, the Bernoulli numbers B n are 
defined as the coefficients of A -1 : 

n 00 

D _ x = ^TB k D k /k\, 

e 1 k—0 


we have 

Po^ix) = (x)^ = log(x + 1) + B\/(x + 1) — B2/2(x + l) 2 + Bz/?>(x + l) 3 -. 

Note that (x — 1)q^ = ipix) is the Gauss psi function. In this domain, we verify 
that 

ip(x) = r'(x)/T(x) 

as follows. Begin with the fundamental identity for the Gamma function, 

T(x + 1) = xT(x). 

11 This information is only preliminary since without the residual theorem (Theorem 0 to 
pinpoint the exact value of p— i(a?), the entire sequence might be off by a factor equal to a 
shift-invariant operator of degree zero a-o + aiD + C 12 D 2 + • • •. 
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Now, take the derivative of both sides. 


r'0 +1) = T(x) + xT(x). 


Next, divide both sides by T(x + 1) = xT(x). 

r(x + i) _ i r(x) 

r(x + i) x + r(x) 

Finally, note that 

I = A m 

x r(x) 

Thus, r'(x)/r(x) equals (x)q 1 ^ possibly up to a constant. 

For the logarithmic lower factorial of order (1) and positive degree n, we 
simply employ the higher order Bernoulli numbers 


Pn \ x ) = ( x )n ) = 

k= 1 


^n-k( x + 1 ) 


where 


D 


=D 


- 1 


OO 


= > B kn D«/k\. 


k -0 


Finally, to compute the logarithmic lower factorials of order a ^ (1), merely 
use the projection maps. 

Pn 0) = ( x )n = sk 'Pa,(l){ x )n ) ■ 


To carry the analogy to the end, we note that by eq. [| and part 2 of Theo¬ 
rem |lo|, we still have 

OO 

E a = J2(a)nA n /n\. 

n —0 

From this, we have the logarithmic version of VanderMonde’s identity 


(x + a)“ = 

k—0 


n 

k 


(a) fc (x)“_ fc 


for all integers n. Similarly, we have the logarithmic analog of the Newton’s 
formula 

OO 

f( x )= ( An f( X ))<x( X )n/ nl 

a t^(O) n=—oo 
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We note that this formula constitutes a substantial extension of Noerton’s for¬ 
mula. For example, if fix) = 1/x, then we have for n negative^] 


Thus, 


A n — 

X 


(—n — 1)! 

= ( Fl w-bi),„, 

(—n — 1)! 

= {x)^} l _ 1 /{—n — 1 )! 

A X=— 1 

= (—1)(—2) • ■ • (n)/(—n — 1)! 

= (—l) ra+1 ?r. 


-1 


i = £ 


n =—oo 


N ! 


= E 


n^o (^ + 1 )(^ + 2 )-"(a: + n) 


(15) 


3.5 Logarithmic Sequences of Binomial Type 

As we did with polynomials, we too can now extend the interplay between the 
derivative D and the forward difference operator into a more general framework 
involving all delta operators]^] Instead of sequences of polynomials p n (x ), we 
have logarithmic sequences p“(x) indexed by an integer n and a vector with 
finite support of integers a. Again we have a requirement that p“(x) be of 
degree n , and also a new requirement that it be of order a ; that is to say, 

n 

Pni X ) = E CfeA fc^) 

— OO 


with c n ^ 0. In particular, pn\x) is a sequence of polynomials: one of each 
degree with ( x ) = 0 for n negative. More generally, the (x) each belong 

to C a and in fact form a basis for it. 

12 Obviously, for n nonnegative, ^ A n (^))^ = 0. 

13 As before, a delta operator is a delta series in the derivative; that is a shift-invariant 
operator whose kernel is the set of constants—in this case R. 
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The only two examples of logarithmic sequences which we have seen so far 
are the harmonic logarithm A“(ir) and the logarithmic lower factorial (a;)“. 
However, there are clearly many more; every sequence of polynomials can be 
extended in many ways into a logarithmic sequence. However, we are particu¬ 
larly interested in sequences of polynomials of binomial type. We seek a natural 
definition of a logarithmic sequence of binomial type such that every sequence 
of polynomials has one unique such extension. The following theorem is exactly 
what we are seeking. 


Theorem 11 (Logarithmic Sequences of Binomial Type) Let p n (x) be a 
sequence of polynomials of binomial type. Then there exists a unique logarithmic 
sequence Pn(x) called the logarithmic sequence of binomial type which 

extends p n {x) 

Pn \ x ) — Vn{x) for all integers n, (16) 

obeys an analog of the binomial theorem 


Pn( x + a ) = Y 
k =0 


n 

k 


Pk(a)p°- k (x) 


(17) 


for all integers n and constants a, and is invariant under the skip operator 
skip a p for /3 ^ (0); that is to say, 

ski P a pPZ(x) = <W“( a 0- (18) 


Proof: Existence Let /(D) be the delta operator associated with p n {x) and 
let g( D) be its compositional inverse.^ Now define p“(x) to be what is known 
commonly as the conjugate logarithmic sequence for <?(D). 


:w = E 


k=—c 


(<?( D) fc A“(x)> c 

M! 


■A?(z). 


After briefly noting that Pn{x) is a well defined logarithmic sequence invariant 
under the skip operators, it will now suffice to prove it obeys eqs. |Ig| and 0- 
This will follow immediately from the following two lemmas. 


Lemma 11.1 (Associated Logarithmic Sequences) Let p^(x) be the con¬ 
jugate logarithmic sequence for the delta operator g(D) whose compositional in¬ 
verse is f{D). Then Pn{x) is the associated logarithmic sequence for f(D) in 
the sense that 

(Pn(x))a = S n0 (19) 

f( D )Pn( X ) = NK-lW- ( 2 °) 

14 They are compositional inverses in the sense of composition of formal power series in the 
derivative. That is, f(g(D)) = D = 
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Proof: Eq. [ 11 ] follows directly from the definition of a conjugate logarithmic 
sequence and the augmentation, 

(pS(*))„ = (^(*))/L 0 l! = 5 nO . 


For eq. ^0[ we must first write /(D) in terms of D n : /(D) = J2kLo ak ^ k - 


Then we proceed as follow 

/( D K 0) = 


(g{ D) fc A“(z)) 

E "n, Ag-jfr) 


k,j> 0 


L fc - j'1 ! 


^ (a,-s(D) fc+ -'A2(a;)) 

E „„ /a Agw 


k,j>0 




E -iTTT-^(z) 


k—0 


m 


“ (Dg(D) fe A“(x)) 

E ,,.i, 


fc=0 


Lfcl! 


“ (g(D) fe A“ Hz)) 

NE—— ,E /a AgW 


k=0 

M Pn-lO^.D 


Lfcl! 


Note that in the case a = (0), pi°^(0) = S„o and f(D)pn\x) = np^li(x), so 
that Pn\x) is associated to the delta operator /(D). Hence, we have eq. |l6|. 


Lemma 11.2 Let Pn(x) be the associated logarithmic sequence for the delta 
operator f(D). Then for all shift-invariant operators h(D) we have 


OO 

H D )Pn{x) = 

k =—oo 


n 

k 


{h{D)Pk(x)) a p*_ k (x). 


Proof: It suffices to prove the lemma for h( D) = /(D)- 7 since the powers 
of the delta operator /(D) form a basis for the space of shift-invariant oper¬ 
ators. However, by the definition of an associated logarithmic sequence, we 
immediately (/(D)- , p^(x)) = 5jk [/]!• Thus, 


/(D) J p“ (x) = 


L n ~|! 


L n - j 1 ! 


E 


k=—oo L 


n 

k 


(f( D yPk(x)) a p“_ k {x).n 
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Finally, eq. [Tl] follows from the lemma with ft.(D) = E a . 

Uniqueness Now it remains only to show the uniquness of logarithmic 
sequences of binomial type. We will do so in two lemmas. The first shows that all 
logarithmic sequences of binomial type are associated logarithmic sequences, and 
the other shows that every delta operator has a unique associated logarithmic 
sequence. These two lemmas will complete the proof of Theorem [H]. 


Lemma 11.3 Let p^(x) be a logarithmic sequence of binomial type; that is, let 
it obey eq. 0 Then there is a unique delta operator HD) such that (x) is 
associated with f{D). That is, Pn( x ) obeys eqs. [7J and\2(\. 


Proof: Eq. ^9| immediately follows from Proposition [I]. 
linear operator Q by 


QPn 0) = \H\Pn- iW- 


Now, we define a 


Eq. ^ is equivalent to the proposition that Q is a delta operator. Since Q 
obviously commutes with skip a/3 for (3 ^ (0), and the kernel of Q is clearly the 
set of real constants R, it remains only to show that Q commutes with the shift 
operator E a . 


QE a p“( x ) = E 

k =0 
oo 

= E 


k=0 L 

oo 


= E 


n 
k 

n 
k 

n — 1 


Pk(a)QPn-k( X ) 

[n — k] p k {a)p a n 

M Pk{a)Pn-k-l( X ) 


k =0 

= E?Qi%{x).U 


Lemma 11.4 Let f(D) be a delta operator. Then there is exactly one logarith¬ 
mic sequence pH x ) associated with f(D). 


Proof: By Taylor’s theorem (Theorem ||), Pn(x) is determined by its aug¬ 
mentations ak = . However, since /(D) fc is a basis for the space of 

shift-invariant operators, the augmentations ak are determined by the augmen¬ 
tations bk = (/(D) fe p“ (x) ^ ^ which we know are equal to kV n / ;; . C QED 

We also derive—in the same manner as Lemma |11.2| —generalizations of the 
logarithmic Taylor’s theorem and expansion theorem. 
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Theorem 12 (Logarithmic Taylor’s Theorem) Let p„(x) be the logarith¬ 
mic sequence associated with the delta operator f(D). Then any logarithmic 
series p{x) G C + can be expanded 


N a 


p( x ) = J2 J2 

a^(0) n= — oo 


(f{D) n p(x)) c 

IN! 


NON 


Theorem 13 (Logarithmic Expansion Theorem) Let Pn(x) be the loga¬ 
rithmic sequence associated with the delta operator f(D), and let a ^ (0). Then 
any shift-invariant operator g{D) can be expanded 


sO) 


(9{D)Pn(x))a 

— — r ^ ^ 


3.6 Explicit Formulas 

We note that if any term of the logarithmic sequence of binomial type pn\x) 
associated with the delta operator /(D) is known a priori, then all of the other 
terms can be computed in terms of it by appropriate use of the various powers 
of /(D) and the projections skip^-^. In particular, we only need to determine 

the term p^\ ( x ) which is called the residual term , and then 

N 0) = ski Pa,(i)/( D ) _ " _1 P-i(a:)/ INI" 


Thus, the following result is extemely fundamental, 


Theorem 14 (Residual Term) Letp^(x) be the logarithmic sequence associ¬ 
ated with the delta operator f(D) : and let a ^ (0). Therf^\ 

P ( -l(x) = — . 

X 

Proof: We will actually prove a seemingly stronger result 

_ Pn( x ) = N !/ / (D)/(D)- 1 ~ ri ^- 1 ( a ')- (21) 

15 The operator /'(D) called the Pincherle derivative of /(D) represents the operator g(D) 
where g{t) = . Thus, if /(D) = D n , then /'(D) = nD n_1 . Another characterization of 

the Pincherle derivative will be given later. (Theorem g) 
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Let g“(a;) denote the right hand side of eq. gl]. Clearly, g“(x) is a logarithmic 
sequence, so by Lemma 11.4, it will suffice to show that g“(x) is associated 
with the delta operator /(D). Of the two defining properties, eq. E is trivial to 
demonstrate. 


/( D k“(z) = In] [n — 1] !/ , (D)/(D) n X1 1 (x) = [n] q%_ i(x). 

To prove eq. [h], we first consider the case n / 0. 

<?£(*)>» = (L n l !/ , (D)/(D) _1_n A“ 1 (x)) a 
= Xl^x)) 

—n \ /a 

However, ^(/(D) - "/ A“ 1 (x)^ is given by the coefficient of D _1 in (/(D) _n ) / ; 
however, this coefficient is always zeroj^] Thus, ( p%(x)) a = 0 for n ^ 0. 

Whereas for n = 0, we have 

(Po(*)) a = (/'(D)/(D)- 1 A/ 1 ( a :)) Q . 

However, /'(D)/(D) _ 1 is equal to D _1 plus terms of higher degree, so 


(Po(*))a = !- D 


Theorem E allows us to derive a series of formulas—as explicit as possible— 
for the computation of sequences of polynomials of binomial type or logarithmic 
sequence of binomial type associated with a delta operator /(D). 


However, we require one last definition. We define on the logarithmic algebra 
the Roman s/ni/f^ as the linear operator a such that for all a 


( A“ +1 (x) for n ^ — 1, and 
\ 0 for n = —1. 


The Roman shift is the logarithmic generalization of multiplication by x. In 
particular, note that for n/— 1, 


We observe the following remarkable facts about the Roman shift. 

16 Thi s theorem will be seen to be equivalent to the Lagrange inversion formula (Theo¬ 
rem |l 7[ l. The observation here that the residue of the derivative of a formal power series is 
zero is the key step in all known proofs of the Lagrange inversion formula. 

17 Named after Steve Roman. Sometimes also called the Standard Roman Shift or the Roman 
Shift for D or the Roman Shift for (x). 
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Proposition 2 (Pincherle Derivative) For all logarithmic series p(x) £ C, 


(Da — aD)p(x) = p(x); (22) 

in other words, 

Da — aD = I. 

More generally, for all shift-invariant operators f(D ) , we have the operator iden¬ 
tity 

f(D)a-af(D) = f'(D). (23) 


Proof: It will suffice to consider the case /(D) = D k and A“(x). When k = 0, 
eq. G3 is a triviality. Now, consider the case in which n is neither k — 1 nor —1: 


(DV — aD k )\“(x) = D fc A“ +1 (x) — 


W! 


\n — k ]! 
[n + 1]! L n l ■ 


crA"_ fe (x) 


A° 


\n - k + 1]! [n - k ]! J n ~ k+1 

N ! 


(x) 


= (L n + l"l - Ln-fc + il) fc + i ^ Ag_ fc+1 (a;) 


= kD k ~ l \ a n (x). 

Next, consider the case in which n = — 1 but k ^ 0: 

(D fc cr — aD fc )A“ 1 (x) = -od^A^x) 


1 


L-fc-1]! 


= ~l-k] 


1 


L-fcl! 


°'A“ fc _ 1 (x) 


= A:D fc_1 A/ 1 (x). 

Finally, consider the case in which n = k — 1: 


(D fc cr — crD fe )A“(x) = D fe A“ +1 (x) — 


L n l- 


[n — k ]! 

= In+il^oO) 

= [n + l]D n \«(x) 

= fcD fe - 1 A“(x).D 


a\f 1 (x) 


Thus, we see that the Roman shift gives a non-trivial extension of the classi¬ 
cal commutation relation Dx — xD = / of quantum mechanics in the logarithmic 
algebra. In effect, eq. E3 is really quite incredible, for Van Neuman has shown 


28 










[ pl| that for polynomials^ the only operators ( 9 , </>) obeying the commutation 
relation. 

6(f) — 4>0 = I 

are 9 = D and <f> = x. However, by eq. ^2], we see that for logarithmic series, 
9 = D and <f> = a is yet another solution. 


Theorem 15 (Rorigues) Let p{x ) be a finite linear combination of the har¬ 
monic logarithms. Then for all real numbers a, we have the following formal 
identity 

(exp(— aa)Deyip(aa))p(x) = (D — al)p(x).0 


We can now state our main result. 

Theorem 16 Let f(D) = Dg{D) be the delta operator associated both with 
the logarithmic sequence p“(x ), and with the sequence of polynomials p n (x) = 
Pn\x). Then we have, 

1. (Transfer Formula) For all n, 

p“(x) = f {D)g{D)~ n ~ 1 \ c f l (x), 

2. (Shift Formula) For n neither 0 nor —1, 

Pn{x) = ag{D)~ n A“_!(x), 


and 

3. (Recurrence Formula) For n ^ — 1, 

Vu{x) = af'(D)- 1 p%_ 1 (x). 

In particular, for the polynomial sequence p n {x), we have 

1. (Transfer Formula) 

Pnfx) = f\D)g{D)- n - 1 x n , 

2. (Shift Formula) 

Pn(x) = X3(R) _n X n_1 , 

and 

18 Since then others have shown that this result holds for much more general sorts of functions 
although obviously not for logarithmic series. 
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3. (Recurrence Formula) 


P„(x) = xf'(D) i ( z ). 

Proof: (1) From eq. |2l], we have for a ^ (0) 

p a n {x) = Lnl!/'(D) 5 (D)” +1 D- 1 -A“ 1 (x) = /'(D) 5 (D)" +1 A“(x). 
For a — (0), we derive the result from the case a = (1) as follows]^] 
p n (x) = skip (0)i(1) pW( 2 ;) 

= s ki P(0) , (1) / , (D)sf(D)" +1 A^ 1 ^(a;) 

= / , (D)g(D)" +1 skip (0)(1) A^ 1 -* (x) 

= f'(D)g(D) n+1 x n . 

(2) For this part, we need the following routine calculation 

/'(D)g(Dr +1 = 5(D)" - ±( ff (D)")'D. 

n 

Thus, by part (1), we have 

KOr) = g(D) n XZ(x) (ff(D) n )'A“_ 1 (a;). 

Finally, by Proposition ||, we have 

p«(x) = 9 (D)»A“( a: )- S (D)VA“_ 1 ( 3: ) + ( T S (D)’>A“_ 1 ( a: ) 

= </(D) n A"(x) - g(D) n A“(a;) + ag(D) n X^_ 1 (x) 

= ag(D) n X^_ 1 (x). 


(3) This follows immediately from the substitute of the Transfer Formula 


A“_ 1 (a:) = /'(D)- 1 ff (D)>«_ 1 ( a ;) 


in the Shift Formula 

Pn( x ) = <? 9 (V)~ n X%_ 1 {x).n 

Actually, Theorem |l6| can be generalized to cover the case in which the role 
of the harmonic logarithm A“(a;) can be played by any logarithmic sequence of 
binomial type (x ). 

19 This is then the best possible example of a purely logarithmic proof of a result involv- 
ing polynomials. It is derived from a formula (eq. p]J) which does not even make sense for 
polynomials. 
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Corollary 2 Let f(D) and g{D) be delta operators associated with the the log¬ 
arithmic sequences Pn(x) and g“(a;), and the sequences of polynomials p n (x) = 
Pn\x) and q n {x) = qn\x). Then 


1. We have 


f'(D)q(D ) n+1 

v a ( x ) = J Vm ’ _ a a (x) 

Pn[ ’ g'(D)f(D) n +' qn[ ’ 


for all a and n, and in particular for a = (0) we have 

, s f\D)g{DT+ 1 , s 

Pn{ ) ~ g'(D)f{Dr+^ n[ - X) ' 

2. And for the polynomial sequences onlyj^\ we have 

p n (x) = x (g(D) n f{D)~ n [x _1 g n (ar)]) . 


Proof: Part one is immediate from the Transfer Formula. Part two results 
from the Shift Formula since the Roman shift of a polynomial is multiplication 
by x.U 


3.7 Lagrange Inversion 

The explicit formulas^] of the preceeding section lead to an elegant proof of 
the Lagrange Inversion formula for the compositional inverse /( _1 )(t) of a delta 
series /(f). 


Theorem 17 (Lagrange Inversion) Let f{t) be a delta series, and g{t) be 
any Laurent series of degree d, then 

= £ / g{t)f\t)f{t)-'- k \) t k . (24) 

fc=d ' X1 (i) 

In particular, 


OO / 

/ M) (*r = £ (*"/'(*)/(*) 

k=n ' 


-1 -k 


t k . 


(i) 


In other words, the coefficient of t k in f( 1 \t) n is the coefficient of t n 1 in 


20 The Roman shift is not in general invertible. 
21 In particular, eq. B 
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Proof: It will suffice to demonstrate eq. [24]. By eq. formula, 

Lfcll^D)^))^ = ^(D)/'(D)/(D)- 1 - fe 

where p £ ( x ) is the logarithmic sequence of binomial type associated with the 
delta operator /(D). Now, 

</(D)>£(*)) q = 5 nk [n ]! = (D"A£(a;)> a , 

so in an augmentation of an operator acting on a p k (x), one can substitute 
/(~ 1 )(D) for D in the operator provided one also substitutes X k (x) for p k {x). 
In particular, 


L^l! (s , (D)Pfc 1) (a:)) (i) = L fc l ! \5(/ ( ^(D))^ 1 ^))^ ■ (25) 

However, by the Expansion theorem (Theorem [id]), [k ]! (g(/^ _1 ^(D))p[ 1 ' ) (a;) 


(i) 

is the coefficient of D fc in D)). Putting all this information together, we 

have eq. 24.□ 


4 Examples 

4.1 The Logarithmic Lower Factorial Sequence 


We now return to our example of a logarithmic sequence of binomial type (Sec¬ 
tion 3.4). The logarithmic lower factorial (x)“ is the logarithmic sequence of 
binomial type associated with the delta operator /(D) = A = — I. 


By Theorem |b|, to determine the residual series, we must first compute the 
Pincherle derivative 

/'(D) = A' = e D = E 1 . 

Thus, the residual series is 


( a;) W=S 1 - = 

x 


1 

X + 1 


This completes our calculations of section 3.4, since all of the other terms (x) 
can be computed via the formula 


a. 

n 


0 ) 


a 

n 


N ! 1 skip a ( i ) D n 1 
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Howver, we should note that in some cases the recursion formula is much more 
useful 

( x )n = a ( x ~ !)n-i- 

In particular, for a = (0) and n > 0 or a = (1) and n < 0, this is equivalent to 

0)“ = x ( x ~ l)n-l- 

In these cases, we immediately have 

(x) n = x(x - 1) • • • (x - n + 1) 

for n positive, and 

f x )W =_ \ _ 

y ,n (x + l)(x + 2) ■ ■ ■ (x — n) 

for n negative. 


4.2 The Laguerre Logarithmic Sequence 


Our next example, the Laguerre logarithmic sequence L^(x) and the Laguerre 
polynomials L n {x) = Ln\x), are the logarithmic sequence and polynomials 
sequence of binomial type associated with the Laguerre operator /(D) = W 
defined by eq. |l0] or eq. |ll|. 


We have for all n and a by the Transfer formula, 


K Or) = /'(D) 


(m) 


K(x) 


= -(D -I) n ~ l K(x) 


= E(-D 


n-\-k 


k—0 


- 1\ [n ]! 


k J [n — k ]! 


A “_*(*)• 


In particular, for a = (0) and n positive, we have 


(26) 


n— 1 

L n (x) = J2 oo(—l)" +fc 

k =0 


n — 1 
k 


7 .! 


t ,—k 


n — k\ 


Note that for the Laguerre logarithmic sequence, we not only have 
ski P(o),(i)^'n ) ( a; ) = L n (x) 
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(27) 


as required by eq. |l8[ but we also have 

ski P(i),(0)Mz) = L n\ x ) 

for n positive.^] This is a very special property; it is only true of logarithmic 
sequences associated with delta operators of the form aD/(D — b) for a, b real 
scalars (o/ 0). 


However, 


L a 1] = log a; + - - \ 
x x z 


which is not equal to skip/,) ^Lq(x) = skip^-, ^1 = log a;. 


From eq. |26| and the classical identity (for polynomials) 

e x D n e~ x = (D - I ) n , 


we have the Rodrigues formula for the Laguerre sequence 


L n {x) = —e x D n ~ 1 e~ x x n . 


(28) 


4.3 The Abel Logarithmic Sequence 


We continue our series of examples with the logarithmic generalization (a;; b ) 
of the Abel polynomials A n (x; b) associated with the delta operator /(D) = E a D 
mentioned in Section O- 


By the shift formla, for n neither 0 nor 1, we have 
K{x) = crX^_ 1 (x — nb). 

For a = (0) and n > 1 or a = (1) and n < 0, this immediately gives us 

x(x — nb) n ~ 1 . 

In particular, the residual series is 

A^(a;) = x(x + b)~ 2 . 

22 The logarithmic sequences of the form a n X^(x) (for example the harmonic logarithm) are 
even more special since they obey eq. [27] even for n = 0. Nevertheless, there are no logarithmic 
sequences satisfying eq. ET] for any negative integers n. 
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It remains now only to compute the terms of degree zero and one. The 
series of degree zero is computed via the transfer formula, and turns out be 
much simpler than might have been expected: 

A^x) = (DE b )'E~ b logo: 

= (I + 6D) log a: 

= logx + b/x. 


As an application, we infer from the logarithmic binomial theorem (eq. 0) 

OO 

A£(:r+ a)=^ 


k =0 L 


0 

k 


the remarkable identity 
b 


x + a 


+ log (a: + a) = - + log x + ^ 


A k {a ; b)A^\x) 

(—l) fe+1 a(a — bk) k ~ 1 x 


fc=l 


k{x + bk) k+1 


For example, we can substitute a = 1, b = 2, and x = 5, and compute the first 
twelve terms of both sides yielding equality to seven decimal places. 

In general, the transfer formula gives us 

AZ(x) = E- nb (I + bD)\«(x) 

= ^“(a; - n&) + & L«l K-i( x ~ nb )- 

For example, the term of degree one is given by 

A^ (x) = x log(x — a) + a — x. 


Finally, consider the Taylor’s theorem for the logarithimc Abel sequence 
(Theorem |l2| ). Any p(x) £ C + can be expanded in terms of the logarithmic 
Abel sequence as follows 

N ° „<* 

= e E 

or^( 0) n=-oo L 1 

where a“ = (E nb D n p(x)) a . In particular, for p(x) = log a;, a“ = 0 for a ^ (1) 
or n positive, but for n nonpositive and a = (1) we have 

c a = (E nb D n logx) {1) 

= (E nb \^(x)/(-ny) w 
= (E nb \W(x)/(-ny)^ 

= {{x + nb)~ n ) {0) /{-n)\ 

= ( nb)~ n /(-n)\. 
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Thus, 


logo; 

b 

x 

b_ 

x 2 


log x h— b i. 

x ' n! —ral! 

n=l L 1 




“ n 

n=l 

A ( nb) n . 

2l - —(x + nb) 

“ n 

n =1 


- 2/2 + nb) 


—n —2 


(29) 


Note that eq. ^9] is the correct version of the calculations in ]l2], p. 105]. We 
thank Richard Askey for bringing our error to our attention. 


4.4 The Logarithmic Ramey Sequences 


Let and p n (x) = Pn\x) be the logarithmic and polynomial sequences of 
binomial type associated with a delta operator /(D). The sequences p"(x-,b) 
and p n (x; b) associated with the delta operator E b f( D) are called the logarithmic 
and polynomial Ramey sequences^ of p n {x) and Pn(x). See Ramey’s paper on 
polynomial Ramey sequences for /(D). For example, if /(D ) = D, the Ramey 
sequences are the Abel sequences dealt with in sections |l.4| and LE. 


The following proposition is very useful in the calculation of Ramey se¬ 
quences. 


Proposition 3 Letp^fx) andp n {x) be the Ramey sequences for the logarithmic 
and polynomial sequences of binomial type q^ix) and q n (x) = q„(x) for the 
delta operator f(D). 


1. Then for n neither zero nor one, the explicit relationship between p n [x) 
and q n (x) is given by 


Pn{x) 


xq n {x — nb) 
x — nb 


2. For all n, we have 

Pn(x) = q%{x - nb) + [n] bf {D^q^^x - nb). 


23 Also called the Abelization. 
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3. For n negative and a = (1), we have 

Pn\x) = q { n\x - nb) + | n] bx~ l q^(x). 

Proof: (1) Immediate from part two of Corollary [|. 

(2) Immediate from part one of Corollary ||. 

(3) Since the Roman shift is multiplication by x in this context, we have by 
the Recurrence Formula 


x 1 Qn( x ) = f'(B)q^_ 1 (x). 

The conclusion now immediately follows from part two.D 

For example, let us now compute the polynomial Rainey sequence for the 
lower factorial. This is the polynomial sequences G n (x) commonly called the 
Gould polynomials associated with the delta operator DA. We have 

x(x-nb) n 

G n (x) = --— = x(x - nb - l) n —i- 

x — nb 


5 Connection Constants 

Let Pn{x) and < 7 “(a;) be logarithmic sequences with 

n 

Qn( X ) = C nkK( x )- ( 30 ) 

k =—00 

Then the sequence 

n 

Xn(x)= C nkPk(x) 

k=— 00 

is called the umbral composition of the logarithmic sequences ( x) and ( x ). 

To denote umbral composition, we write 

C{x) =g“( p). 

Restricting out attention to polynomials, we see the if p n (x) and q n {x) with 

n 

Q.nip^) — ^ ^ Cnkfi 7 
k—0 
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then the umbral composition is uniquely defined, and is given by 

n 

q n { p) = y,c nkPn {x). 

k =0 

When (x) and < 7 “ (a;) are of binomial type, then umbral composition has 
a remarkable property. 

Theorem 18 Let Pn(x) and qn(x) be the logarithmic sequences of binomial 
type associated with the delta operators f(D) and g(D) respectively. Then their 
umbral composition r™{x) =p“(q) is the logarithmic sequence of binomial type 
associated with the delta operator f(g(D)). 

Proof: We must show the r“(a;) obeys eqs. ^9] and |o|. Let the coefficients of 
Pn(x) and /(D) be given by 

n 

Pn(x) = y C nk \l(x) 

k— — oo 
00 

/(D) = 5>,ry 

j =1 

We write c nk since by eq. |l 8 | the coefficients do not depend on a. 

To prove eq. |l^, we proceed as follows: 

(C(q))« = 


( y c nk q%(x) 

\ k= — oo / 

n 

y c nk {q%{x)) c 


k =—00 
n 

^ ^ Cnk^Ok 
k =—00 
^nO 

( Pn( x )) a 
dnO • 


Next, it remains to prove eq. ^0|. 

00 n 

f(g( D))r“(x) = y y ajC n kg(Dyq%(x) 

j =1 k=—oo 

°° n I I 

= E E a jCn k 1 \ v q k -j{x). 

j=l k=-oo 
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However, 


L n l Pn-l{ X ) = f(D)Pn(x) 


oo n 


= J2 a j c nkV J \k(x) 

j =1 k =—oo 

Lfcl! 


oo n 


= a i Cnk 


K-j 0 * 0 , 


j—1 k=—oo 




so that 


/(s(°)KW = M C~i( x )- D 


By restricting our attention to we have the analogous result for poly¬ 
nomials. 


Corollary 3 Let p n (x) and q n {x) be the polynomial sequences of binomial type 
associated with the delta operators f(D) and g(D) respectively. Then their um- 
bral composition r n (x) = p n ( q) Is the sequence of binomial type associated with 
the delta operator f(g(D)).D 


The preceeding theorem allows us to solve the problem of connection con¬ 
stants. That is, given two logarithmic sequence of binomial type r“ ( x ) and (x) 
associated with the delta operators g(D) and h(D), we are to find constants)^] 
c„.k such that 

n 

r n( X ) =J2 Cnkq kW- 
— OO 

The preceeding theory shows that the constants c n k are given by a logarith¬ 
mic sequence p„{x) as in eq. 1 associated with the delta operator /(D) = 
/i(<j/ -1 )(D)). Of course, the same result holds true for polynomial sequences of 
binomial type. 

In particular, given any logarithmic sequence of binomial type Pn{x) asso¬ 
ciated with the delta operator /(D), there exists a unique inverse logarithmic 
sequence p%( x ) associated with the delta operator /1 _1 )(D). In view of the 
logarithmic Taylor’s theorem (Theorem |j), 

n 

Pn( X ) = 

k=—=infty 

24 Such constants exists since the two sequences are basis, and the coefficients do not depend 
on a by eq. hq. 
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However, by the demonstration of eq. g5|, p“(x) is the conjugate logarithmic 

sequence for the delta operator /(D); that is, 



k =—oo 


Moveover, by further umbral composition 



k =—oo 


Obviously, similar results hold for polynomials. 


6 Examples 

6.1 From the Logarithmic Lower Factorial to the Har¬ 
monic Logarithm 

To compute the connection constants expressing the harmonic logarithms A“ (x) 
in terms of the logarithmic lower factorials (x)“ via the above techniques, we 
must first calculate the logarithmic generalization </>“ (x) of the exponential poly¬ 
nomials 4> n {x) associated with the delta operator A^ 1 ) = log(I + D) and which 
have been studied by Touchard Q among others. 

As noted above, this is the logarithmic conjugate sequence for forward dif¬ 
ference operator A. That is, 



k =—oo 


In particular, for a = (0) 



where S(n,k) denote the Stirling numbers of the second kind. 
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To calculate the residual series, we as usual employ Theorem [l4|, 

= log(/ + D/a ;" 1 

OO 

= l) fc D fe :c _1 

k =0 

OO 

= J2 k] / xk+1 

k =0 

which is equivalent via umbral composition to eq. [u|. 

Other terms of the logarithmic exponential sequence are most easily com¬ 
puted via the recurrence formula. 

CO) = Cr(log(/ + D )') _ 1 C-lO) 

= a(I + D)C-i(x). 

For a = (0) and n positive, or a = (1) and n negative, we thus have by eq. ^ 8 ] 

CO) = ze _a: De x C_iO) 

= e-*(:rD)e*C-iO). 

In particular, for polynomials, we have by induction 

cj> n {x) = e- x (x D) n e x . 


6.2 From the Lower Factorials to the Upper Factorials 

The logarithmic upper factorial (a:)“ is the logarithmic sequences associated with 
the backward difference operator V = I — E~ D . In effect, it is the important 
special case of the Gould sequence with b = — 1. Thus, for a = (0) and n 
positive, we have the sequence of polynomials 

0 )n = 0)n = x(x + 1) • • • 0 + n - 1), 

and for a = (1) and n negative, we have 

(C 1 ) =_ 1 _ 

V ,n {x — l)(x — 2) -■- {x+ n)‘ 


To express (a:)“ in terms of 0)0 we m ust first compute the logarithmic 
sequence of binomial type p“0) associated with the delta operator /(D) = 
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A(V(-^). 


/(D) = A(V (-1 )) 

= exp(— log (/ — D)) — I 


D 

I -D 
= -W 

where W is the Weierstrass operator. Hence, the logarithmic sequence p“(ai) is 
essential the logarithmic Laguerre sequence. 

rtW = (-I)’KM = D- 1 )^" t 

Thus, 

In particular, for a = (0) and n positive, 

x(a; — 1) • • • (x — n + 1) = E( _1 ) fc ( n ) — —n x ( x + 1) • • • {x + n - k - 1) 

a fact which was not at all obvious a priori. 


6.3 From the Abel Series to the Harmonic Logarithms 

To compute the harmonic logarithms A “(a?) in terms of the logarithmic Abel 
sequence A“hr; b) it will suffice to compute the logarithmic inverse Abel sequence 
b) = A n (x', b) associated with the composition inverse of D E b . 

Unfortunately, there is no closed form expression for the compositional in¬ 
verse of Dexp(6D). Nevertheless, we can compute p^ix) using the fact that it 
is the conjugate logarithmic sequence for D E b . 

n (D k E bk \%(x)) 

<(*)= E A -Fit— 

Now, (D k E bk X^(x)^ is equal to ( bk) n ~ k \n\\/\n — fc]!. Thus, 

<M> = ± 
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And therefore 


Ki x ) = 


m 


n—k 


k=—oo 


N! 


-Ak{x-,b), 


or for a = (0) and n positive or a = (1) and n negative, 

" ( bk) n ~ k . 

x = 2^ x(x-nb) n 


k— — c 


N! 
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